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ABSTRACT 

We study the 0509-67.5 supernova remnant in the Large Magellanic Cloud with the VLT/F0RS2 spectrograph. 
We detect a broad component in the Ha emission with a FWHM of 2680 ± 70 kms^^ and 3900 ± 800 kms^^ for 
the southwest (SW) and northeast (NE) shocks respectively. For the SW, the proton temperature appears to be too 
low for the shock velocity, which we attribute to a cosmic-ray pressure behind the shock front of at least 20% of the 
total pressure. For the NE, the post-shock proton temperature and the shock velocity are compatible, only if the 
plasma behind the shock front has a degree of thermal equilibrium of over 20%, which is at odds with current models 
for temperature equilibration behind fast shocks, which do not accelerate cosmic rays. If we assume the electron 
temperature to be less than 10% of the proton temperature, we find a post-shock cosmic-ray pressure of at least 7%. 
Subject headings: ISM: individual objects (SNR 0509-67.5) — ISM: supernova remnants — radiation 
mechanisms: thermal — acceleration of particles — shock waves 
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1. INTRODUCTION 

Supernova remnants (SNRs) are generally thought to 
be the dominant sources of Galactic cosmic rays. For this 
to be true, SNRs need to transfer about 10% of their ini- 
tial kinetic energy into cosmic rays. An open question 
is whether SNR shocks can reach these acceleration ef- 
ficiencies averaged over their life time; recent TeV and 
GeV 7- ray observations give promising but ambiguous 



results (Ellison et al.||2010 Abdo et al. 2010) 



Although the process of efficiently acceleratin g parti- 
cles is well understood (Malkov & Drury 2001), obser- 
vatio nal verifications for efficient acc eleration are scarce 
Warren et al.|2005l|Lee et al.|2007; Vink et al. 2006; 



Helder et al. 2009). These observations are essential for 
characterizing the efficiency of the acceleration. As the 
particles move ahead (i.e. upstream) of the shock, they 
form a so-called shock precursor, which compresses and 
pre-heats the upstream medium. This effectively lowers 
the Mach number of the main shock, there with lowering 
the temperature behind the s hock front (Ellison et al. 
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O'C. Drury et al. 2009). This effect would ide 



e characterized by the post-shock proton tempera- 



ture (Tg p), as Tg p is close to the plasma temperature, 
whereas the po st-shock electron temper ature Tg^ e rnight 



be lower (e.g. Ghavamian et al.||2007b ). The latter im 



plies that the electrons might constitute only a minor 
part of the thermal pressure behind the shock front. 

For some SNRs, Tg, p can be determined from hy- 
drogen line emission at the shock fronts; these are so- 
called Balmer-dominated shocks. The hydrogen lines 
of a Balmer-dominated shock consist of two superim- 
posed components; the narrow component is emitted by 
neutral hydrogen after entering the shock front and the 
broad component by hot protons after undergoing charge 
exchange with incoming neutral hydrogen atoms. The 
width of the broad component re flects the proton tern 
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perature be hind the shock front (Chevalier et al 
' Henil[2009|) . 

In this Letter, we report on our study of SNR 0509 

Electronic address: |e.a.helder@astro-uu.nl| 



67.5, as this is a likely source of efficient cosmic-ray ac- 
celeration; recent studies show that the re mnant resulted 
from a highly energetic Type la explosion ( Hughes et al. 



1995| [Badenes et al.|[2008l [Rest et al]|200 8). Its X-ray 



spect rum provides some evid ence tor non-thermal emis 



sion ( [Warren fc Hu ghes 2004[) 
ical models show 



In addition, hydrodynam- 
that "the pres ence of highly energ etic 



particles is likely not negligible ( Kosenko et al.[[2008 ). 
SNR 0509-67.5 was discovered to have Balmer 



dominated shocks by Tuohy et al. (1982), although the 
flux in the broad compone nt was too low to be detected. 

did not 



A subsequent attempt by Smith et al 
reveal any broad component either 



( |1991[ ) 
Tlie first detection 



of a broad compon ent in the hydrogen line emission of 



this remnant was by [Ghavamian et al. (2007a), who mea- 
sured the width of the broad component of the Ly/3 line 
to be 3700 ±400 kms~^, corresponding to a shock veloc- 
ity {vs) of 5200-6300 kms^^. However, as the spectrum 
was taken from the entire remnant, it remains uncertain 
whether the obtained lincwidth is broadened by the bulk 
motion of the plasma as well. 

We report the detection of a broad component in the 
Ha line emission of SNR 0509-67.5 at two locations of 
the shock front. We combine this with a shock veloc- 



ity based on X-ray observations (Kosenko et al. 2008), 
to determine the fraction of the post-shock pressure con- 
tributed by cosmic rays. 

2. OPTICAL DATA AND RESULTS 

We observed SNR 0509-67.5 for 10932 seconds (4 
blocks of 2733 seconds) with F0RS2, the low dispersion 



spectrograph of ESO's VLT ( [Appenzeller et al.||1998 ) 
The observations were made on October 15 and 20 and 
November 10 and 11, 2009. By centering the slit at a 
bright star at a = 05 : 09 : 28.793, S = -67 : 31 : 30.83 
(J2000), with a position angle of 42° (Figure [I]), we ob- 
tained spectra of both the SW and NE rims with a sin- 
gle pointing. The width of the slit was 1.6" which, in 
combination with the 600RI grism and the 2x2 binned 
readout, corresponds to a resolution of ^ 485 kms~^ at 
Ha (6563 A). This resolution prevents us from resolving 
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Fig. 1.— Left: Spectrum and best fit of the NE shock of SNR 0509-67.5. Middle: Archival HST/ACS observation of SNR 0509-67.5, 
as obtained at October 28, 2006, with an exposure time of 4620 seconds. The position of the longslit of the VLT/FORS2 instrument is 
indicated. The star on which the sht is centered is indicated with a circle. North is up, east is to the left. Right: Spectrum and best fit 
of the SW shock. The lower panels show the residuals divided by the errors. The insets show a zoom-in of the broad component, with the 
best fit narrow and broad component overplotted. 



the narrow component of the Ha hne, but it increases 
the signal-to-noise sufhciently to detect and resolve the 
broad component. The data were corrected for bias, flat- 
fielded, and the skylines and cosmic rays were removed. 
The wavelength calibration was done by fitting a fourth 
order polynomial to the spectral lines of He, HgCd, Ar 
and Ne lamps, obtained during daytime. In addition, we 
checked our calibration against the position of 3 skvlines 
(at 5577.3, 6 864.0 and 7571.7 A respectively, 



et al. 1996). The absolute wavelength calibration has 



Osterbrock 



a systematic uncertainty of 0.5 pixel ( 37 kms~^) and 
the errors on the datapoints were determined by trac- 
ing the error propagation through the reduction steps, 
starting from the raw data. For each observation block, 
obvious outliers, probably caused by erroneous skyline 
subtraction, were removed. The resulting spectra were 
fitted with two superimposed Gaussians and an offset, 
convolved with the slit width. The best fitting param- 
eters for both spectra are listed in Table [2J In order 
to determine the significance of the broad component in 
the NE spectrum, we compared the obtained x^ value 
(464) with a fit in which we only fitted a single Gauss 
to the spectrum (525). The difference indicates a sig- 
nificance of 7.8(7. The Xrcd values are 2.03 and 2.56 for 
the SW and NE fits respectively. These high values are 
mainly caused by substructure of the narrow line: a fit 
to the broad component with a single Gauss, excluding 
the central region (between -100 and 800 kms""'^), con- 
volved with the resolution, gives Xrcd '^^ 0.51 for the SW 
and 0.76 for the NE. This substructure might either be 
caused by spatial surface brightness fluctuations of the 
remnant within the slit, or by errors in the skyline sub- 
traction. We determine the la errors on the parameters 
for the broad components by using the errors we deter- 
mined in our data reduction. For the narrow component 
we first increased the errors such that Xred = 1' ^^"^ then 
determined the 1 a errors, using Ax^ = 1. 

3. COMPARING WITH SHOCK VELOCITIES 

To determine the shock velocity, we use the X-ray line 
width of ctrgs — 4900 ± 420 kms~^, as observed by the 
Reflection Grating Spectrometer (RGS) onboard XMM- 



Newton (Kosenko et al. 12008 1. This line width is caused 



by both the thermal and bulk broadening of th e plasma. 



We follow the method of Kosenko et al. (2010) to disen- 



tangle the bulk broadening from the thermal broadening 
to obtain the forw ard shock velocity, a s described below 



As shown in Kosenko et 
related as follows: 



al. (2010), aRGS and Vs are 



CTRGS = (ws/4)y 3r2 



q^2 



Ik' 



(1) 



where Tgh = Wrs/ws^ the ratio ot the reverse shock 
velocity to the forward shock velocity and rbuik = 
^buik. cj/'^^buik, CSM^ the gradient in the plasma bulk ve- 
locity from reverse to forward shock. We obtain rgh — 1 
(whole remna nt) and rsh — 0-5 (NE), from analytical 
models ( [Truelove fc McKee|[l999l ) with Mej = 1-4 Mp , 
E = 1.4 X 10'^^ erg and a n age of 400 years (Badenes et al 



2008| Rest et al 

models with a forwarc 



2008|). Additionally, we constrain the 
shock radius of 15.9 ± 0.8" and 



16.3±0.3" for the entire remnant and the NE respectively 
and a reverse shock radius of 11.4 ± 2.1" and 13.5 ± 0.2" 
for the entire remnant and the NE respectively, based on 
a deprojection of Chandra images (48.9 k s , obs ID:776), 



following the procedure of IKosenko et al 
obtain 



We 



J^bi 



Ik 



0.95 with 



^(|2010) 

numerical 



(Sorokina et al. 2004 Kosenko 2006), 



simulations 
using the above 
parameters. 

We use equation ([I]) to estimate Vs = 6000±300 kms~^ 
(whole remnant) and fs = 6600±400 kms~^ (NE). In the 
remainder of this paper, we use conservative Vg estimates 
of 5000 kms-i for the SW and 6000 kms"! for the NE. 

To check these shock velocities, we estimated the 
expansion with Chandra. Chandra observed SNR 
0509-67.5 for 3 times, with the first observation in May 
2000 (earlier used in this section to determine forward 
and reverse shock radii) and the latter two in May 2007 
(29.5 and 32.7 ks, obs IDs 7 635 and 8554). Following 
the method of Vink (2008b), we find a shock velocity 
of 6700 ±400 kms~ averaged over the azimuth of the 
remnant. Note that a full proper motion study of SNR 
0509-67.5, using both Chandra and Hubble data, is 



underway 
in the SW 



([ Hughes et al.||2010 l 
(I'lg. [I]) suggests 



The dense material 
that the forward shock 





Ha line 


emission 


of SNR 0509-67.5 




Parameter 












SW 


NE 


center narrow 






[kms-i] 






287.0 ±1.4 


286.0 ±1.5 


center broad 






[kms-i] 






-342 ± 28 


459 ± 220 


FWHM broad 






[kms-i] 






2680 ± 70 


3900 ± 800 


total flux 


fio- 


'iSer 


g s~^ cm~ 


arcsec" 


-^1 


8.6 


5.3 


Jb/-fn 












0.29 ±0.01 


0.08 ±0.02 



Approximate flux calibration based on an observation of a photometric standard star (LTT 2415, [Hamuy et al.|[l994^ taken 
November 11, 2010 



might have recently slowed down. However, we do not 
find any support for this scenario from the Chandra 
expansion study, nor from our XMM-Newton/RGS 
study. Moreover, as the RGS study is skewed towards 
the bright SW region, it is unlikely that the SW velocity 
is substantially lower than 6000 kms^^. 



4. INTERPRETATION 

The centroid of the broad component in the SW of 
—342 ± 28 kms^^ indicates the bulk line-of -sight ve- 
locity of the shocked protons. Additionally, the Large 
Magellanic Cl oud (LMC) is movi ng with a 278 kms~^ 

This means that we 
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away from us (Richter et al. 

are observing a part of the shell which is moving toward 
us with 620 kms~^ with respect to the LMC. Figure [l] 
shows a bright, inner shell close to the outer shock. Tak- 
ing into account the seeing of 0.8-1.0", our spectrum is 
probably contaminated by emission from this shell. This 
pa rt of the shell is likely to d ominate the measurement 



by Ghavamian et al. (2007a), as they also measured a 
positive offset for the centroid of the broad component. 
The flux in the broad component with respect to the 
flux in the narrow compone nt (Ih/In ) declines as a func- 



tion of shock vel ocity (Hcng & McCray 2007; van Adels 
berg et ar]|2008 ) . Our low values for Ih/In are therefore 



consistent with a high shock velocity with likely the high- 
est shock velocity in the NE, which has the lowest Ih/In- 
However, as the standard models for interpreting Balmer 
dominated shocks do not include the effects of cosmic-ray 
acceleration, it is not appropriate to use the Ih/In values 
to determine the shock velocities of SNR 0509-67.5. 
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Fig. 2. — FWHM as function of Vb for non-accelerating shocks, for 
different values for Te/Tp, (courtesy to K. Heng and M. van Adels- 
berg). The grey (black) lines indicate this relati on for plasmas, 
that are optically thin (thick) to Ly/3 scattering | |van Adelsberg| 
[eTal. 2008). Overplotted is Vs and the measured FWHM tor both 
We SW and NE. 



4.1. Interpreting the FHMWs 

To determine the cosmic-ray pressure from our mea- 
surements, we need to determine Tg^ p from the FWHM. 
For low temperatures and s hock velocities, kT^^ p = rripa^ 
(JRybicki fc Lightman|l979|, with a = FWHM/ V8 In 2 in 



cm s '-. However, as cross sections for ch arge exchange 



decline for high velocities (e.g. Fig. 1 in Heng & Mc 



Cray 2007), kTs^p > TOpCr^ for higher temperatures and 



shock velocities. Recent studies focused on determining 
the FWHM as function of Vs for n on-acc elerating shocks 



Heng &: McCray|[2607| |vaji Adelsberg et al 
Fig. [2]). Here, we are interested in F 
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Whm 



(e-g- 

resulting m 

as function of Ts_ p for a given Vg instead. We approxi- 
mate this concave function linearly between 0,0 and the 
FWHM and Ts.p expected for a non-accelerating shock 
in thermal equilibrium (Fig. ^and^. In this way, we 
overestimate Tg. p and the corresponding thermal pres- 
sure, leading to a conservative measure of the cosmic-ray 
pressure behind the shock front. 
We note th at the y-axis of Fig. 5 of van [van Adelsberg] 
(2008) is labeled incorrectly (M. van Adelsberg, 
private communication). Instead of showing the 
'broad Ha FWHM', this figure plots the 'broad neutral 
velocity distribution FWHM', which is independent of 
the emission line considered. To model a specific emis- 
sion line such as Ha, one has to convolve the broad neu- 
tral velocity distribution with the relevant atomic cross 
sections. The FWHM-Wg relations we use in this study, 
obtained in electronic form from M. van Adelsberg and 



expected 



measured 




Fig. 3. — Cartoon of a^ as function of Ts^p. 'Expected' are 
the values we expect for a non-accelerating shock (Fig. pi. 'True' 
shows the true Ts_ p and 'measured' is the corresponding a^ we 
would measure. 'Interpolated' shows our estimated Ts^p, based on 
the measured a^ . Note that this method leads to an overcstimation 

0fTa,p. 
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K. Heng, are plotted in Fig. [2] and are based upon the 
same relations used to gene rate Figure 13 and Table 1 of 
van Adelsberg et al.| ([20081 ). 

ink ! ( |2005a| '^and [Helder et aLJ ( [2009[ ), we 
and w, m terms of cosmic-ray pressure be- 



Ibllowmg 
interpret Tg^ p 



hind the shock front and cosmic-ray energy flux leaving 
the system (respectively wcr = PcR/Ptot in which Ptot 
is the total pressure behind the shock front (x-axis in Fig- 
ure |4]) and eosc = F^sc/Ftot in which ftot is the energy 
flux entering the shock; ^pisM^s (Y'^^is of FigureWl). To 
conservatively estimate the post-shock cosmic-ray pres- 
sure, we assume the electrons and ions to be in thermal 
equilibrium. We add wqr and Cesc to the equations of 
conservation of mass, momentum and energy over the 
shock front, which leads to: 

/cTs^p = (1 - it;cR)-(l )iimpvl. (2) 

A /I 

Formally, /i is the number averaged mean particle weight, 
(^ 0.54 for a fully equilibrated and fully ionized plasma 
with LMC abundances), when considering fcTg^p, we can 
treat /x as well as a measure for the temperature equi- 
libration behind the shock front, where fJ- — 1 indicates 
no temperature equilibration, and fi — 0.54 indicates a 
fully equilibrated plasma. In addition, x is the total 
shock compression ratio. For a non-accelerating, adia- 
batic shock, x = 4 and hence kTg.p = jQiinipV^. We 
define /3 = kTs^p/jQfJ.mpVg, to characterize the influence 
of cosmic-ray acceleration on Tg^p. Figure [4[ shows /? in 
the {wcR, eosc)-fi'anie. The 'max'-line indicates the ratio 
of the cosmic-ray pressure and escaping cosmic-ray en- 
ergy for t he most efficiently accelerating shock a ccording 
to theory (|Malkov|1999[|d'C Drury et al.|2009[ ). Hence, 
the hashed region is excluded. 

4.2. The southwest spectrum 

For Vs — 5000 kms^^, we expect to measure at least 
a FWHM for the broad component of 3600 kms~^, cor- 
responding to a Ts,p of 28.7 keV (Fig [2]). This differs 
significantly from the measured 2680 ± 70 kms~^. We 
derive a T^^p of 15.9 ± 0.9 keV, leading to a /3 < 0.58, 
which gives a cosmic-ray pressure behind the shock front 
of at least 20% (Fig. [I]). 

4.3. The northeast spectrum 

For the NE, the measured FWHM indicates a. Vg = 
,/T,,p < 1.0 (Fig. § 



This 



6000 kms-\ for 0.2< T.^e/ 

leads to two possibilities: 

1) The shock does not efficiently accelerate cosmic rays 

and has a Tgc/Tg.p > 0.2, breaking with the earlier re- 

?iOrtcd tr end of Tg e/J's.p oc l/v^ for Vs > 400 kms~^ 
Ghavamian et al. 2007a), as this would give Tc/Tp < 



0.01 for Vs = 6000kms-^ 
2) The shock is far out of thermal equilibrium, as we 
might expect for a fast shock, and is accelerating parti- 
cles. If we assume Tg.c/Ts^p < 0.1, we expect a FWHM 
of >5100 kms~^, whereas we mea sure 3900±800 kms~^. 
Using the approach from section |4.1| and Fig. [4| we ob- 
tain /3 < 0.85 and hence a cosmic-ray pressure behind 
the shock front of > 7% of the total pressure. Note that 
this is a very conservative lower limit: for both Tgp and 
Vs we used conservative approximations. Also, at this 
high shock velocity, the squared line width of the broad 



I.I I I.I 1 1 1 1 J I I.I 1 1. 
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Fig. 4. — Values for /3 in the (uicRi ^csc)-fi'ame. The x-liii6s 
indicate the compression ratio of the plasma behind the shock front. 
The hashed region requires a shock, more efficient than can be 
explained with current theoretical models. The red line indicates 
the lower limit for /3 for the SW. For the NE, /3 is determined 
assuming Ts,c/7s,p < 0.1. 

component does not increase linearly with Tg.p, but fiat- 
tens off. This makes the difference in line width lower 
for a given ATg p, which makes that the line width is 
a less sensitive temperature indicator at high shock ve- 
locities. Future progress can be made with better shock 
velocity estimates of the NE region itself, a higher signal- 
over-noise spectrum of the NE region and with models 
for non-radiative Balmer dominated shocks that include 
the effects of cosmic-ray acceleration. 

A remaining question is whether the magnetic field 
pressure makes a contribution to the post-shock pressure. 
As SNR 0509-67.5 is at a distance of 50 kpc, we can not 
resolve filaments of X-ray synchrotron emission, which 
are often used for determ ining post-shock m agnetic field 
strengths (Vink & Laming 2003). However, Volk et al. 
( 2005 1 showed that a typical value for the magnetic field 
pressure is around 3.5% of the total post-shock pressure. 
Moreover, according to Bell's theory ( |Bell|2004 ) the mag- 
netic field energy density scales as BySiT^ ^VgUc/c, 
with Uc the cosmic ray energy density. This means that 
the magnetic field energy density is expected to be about 
1% of the cosmic-ray energy density. 



5. CONCLUSIONS 

We investigated the cosmic-ray acceleration efficiency 
of the 0509-67.5 SNR in the LMC, by comparing T^^p, 
determined from the Ha-line widths of the SW and NE 
shocks with shock velocities of respectively 5000 kms~^ 
and 6000 kms~"^ for the SW and NE, based on X-ray 
observations. Our study gives the following results: 

- We measured widths of the broad components of 
the Ha fines in SNR 0509-67.5 to be 2680 ± 70 
kms^i for the SW and 3900 ± 800 kms"^ for the 

NE. 

- For the SW, we can only explain the FWHM if we 



Ha line emission of SNR 0509-67.5 



allow for a contribution of > 20% of the post-shock 
pressure by cosmic rays. 

For the NE shock, we have two options: either the 
shock has a Ts^c/Ts^p > 0.2, breaking with the ear- 
lier repo rted trend of Tg, c/^s. p o^ l/^s foi' ^s > 400 



km s~ ^ ( Ghavamian et alj2007a I , or if we assume 
Tg, e/Ts, p < 0.1, the cosmic-ray pressure behind the 
shock front is at least 7% of the total pressure. 

This research, t ogether with our previous study 
( Helder et al. 2009 1 , shows that more than 10% of the 



pressure m young SNRs can be contributed by cosmic 
rays. This is more than the requirement that 10% of 
the available energy needs to be in cosmic rays. On the 
other hand, the cosmic-ray acceleration efficiency may 
decline for older S NRs, as indicated by a recent study of 
the Cygnus Loop ( )Salvesen et al.|2009 ) . So a higher efh- 
ciency at a young age may be needed to have an average 
efficiency of 10% over the whole lifetime of a SNR. 
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